Implications of the CDF tt Forward-Backward Asymmetry for Boosted Top Physics 
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New physics at a high scale A can affect top-related observables at 0(1/ A 2 ) via the interference of 
effective four quark operators with the SM amplitude. The (u , y ll ,'y 5 T a u){t / y >1 ^i 5 T a t) operator modifies 
the large M t j forward-backward asymmetry, and can account for the recent CDF measurement. 
The (u^ tl T a u)(t^ IJ 'T a t) operator modifies the differential cross section, but cannot enhance the 
cross section of ultra-massive boosted jets by more than 60%. The hint for a larger enhancement 
from a recent CDF measurement may not persist future experimental improvements, or may be 
a QCD effect that is not accounted for by leading order and matched Monte Carlo tools or naive 
factorization. If it comes from new physics, it may stem from new light states or an C(l/A 4 ) new 
physics effect. 
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Introduction. The top quark is unique among 
the known elementary fermions in that its coupling to 
the electroweak symmetry breaking sector is not small. 
There is still much to be explored in both the top quark 
sector and the electroweak breaking sector. This situa- 
tion makes the experimental study of top physics inter- 
esting as a probe of new physics, and promising in its 
potential to lead to actual discoveries. The Tevatron ex- 
periments, CDF and DO, are now reaching a stage where 
the precision of their top-related measurements might 
provide first hints to such new physics. 

The CDF Collaboration has recently provided two 
new intriguing measurements. First, a large forward- 
backward tt production asymmetry was observed for 
large invariant mass of the tt system [![: 



A u (M t i > 450 GeV) = +0.475 ± 0.114 , (1) 



to be compared with the Standard Model (SM) predic- 
tion [2- ; 4], Af = +0.09 ± 0.01. 

Second, the CDF Collaboration has recently made 
progress in studying the mass distribution of highly 
boosted jets (p T > 400 GeV for the leading jet) [5{ . This 
study led to an upper bound of 20 fb on the correspond- 
ing boosted top pair cross section, based on naive QCD 
background estimation [f|. The analysis included two 
channels, one involving two massive jets (130 — 210 GeV) 
and another with one massive jet and large missing en- 
ergy. An interesting result found was a significant devia- 
tion from the estimated background in the first channel, 
while no excess was found in the second channel or in the 
combined inclusive search. However, in Ref. 0] it was ar- 
gued that the hadronic channel is more sensitive to the 
presence of boosted tops, and accounting for the excess in 
that channel leads to a tension of less than 1.5 standard 
deviations in the missing Et channel. This observation 
motivates us to consider the possibility that the excess is 
associated with an enhanced boosted tops cross section, 
which might also be linked to Eq. ([1]). 



The estimation of the excess depends on a parameter 
Rmass, described below in Eq. ([3]), which determines the 
QCD background. Assuming that both the statistical 
and systematic uncertainties scale linearly with -R^ssi 
the cross section for ultra-massive boosted jets (not com- 
ing from QCD events) can be written as follows 

cr b = cr th * h (p T >400GeV)- [21-(8.7±3.1)i? m L s ] fb, 

(2) 

where th stands for a hadronically decaying top. The SM 
prediction is crf M = 2.0 ± 0.2 fb @|. 

It is not unlikely that the differences between either 
or both of these measurements and the corresponding 
SM predictions will disappear with improved experimen- 
tal precision, or will be explained by non-trivial QCD 
effects. Yet, either or both of these effects might repre- 
sent hints for new physics. Our approach in this work is 
the following. We interpret the measurement of in 
terms of new physics, checking the consistency of such 
a scenario with other measurements that do not show 
any significant deviation from the SM predictions. Then 
we extract the predictions of such new physics explana- 
tions for ultra-massive boosted jets at the Tevatron, and 
compare to the recent measurement. 

Several works have interpreted the recent CDF mea- 
surement of A^ within specific models of new physics [9j- 
[T5I ]. Similarly, new physics models were invoked flH-tiol 
and model-independent studies were performed [3ll433j 
to explain earlier DO and CDF measurements of the in- 
clusive asymmetry 34j, |35j . 

We do not discuss a specific new physics model, but we 
focus on a large class of models with the following two 
ingredients: 

• The scale of the new physics is well above the scale 
M t i that is relevant to the CDF measurements. 

• The dominant contribution to A£ comes from inter- 
ference between the new physics contribution and 
the SM contribution to top pair production. 
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These assumptions allow us to follow a low energy model 
independent approach, and lead us to particularly clear 
and strong conclusions. Ref. [32j has recently presented a 
comprehensive analysis of top pair production at hadron 
colliders within the same framework. The novelty in our 
work is, first, the incorporation of the measurement of the 
Mjj-dependent A u and, second, the study of the boosted 
jets. 

Boosted jets production. The CDF study [6j fo- 
cused on events with two jets, with a lower bound on 
the transverse momentum (px > 400 GeV) and an up- 
per bound on the pseudorapidity (77 < 0.7) of the leading 
jet. As concerns the jet masses, CDF has defined "light" 
(30 - 50 GeV) and "massive" (130 - 210 GeV) jets. The 
search was divided to four regions. Region A corresponds 
to events with two light jets, regions B and C to one light 
and one massive jet, depending on which is the leading 
jet in terms of pt, and region D corresponds to two mas- 
sive jets. The top pairs should contribute to region D. 
To estimate the QCD contribution to this region, three 
assumptions were invoked: 

1. Events in regions A,B,C come from only QCD; 

2. The actual cross section can be factorized into the 
partonic cross section, which only weakly depends 
on the masses of the final states, and the jet and 
soft functions; 

3. The masses of the leading and sub-leading jets are 
largely uncorrelated. 



Under these assumptions, 



= 1, 



(3) 



where nx is the number of QCD events in region X. As- 
sumption 3 above could turn out to be wrong if there 
is some mechanism in QCD which leads to bias towards 
two massive jets. In Q it was shown that i? maS s is insen- 
sitive to the variation of the relative partonic momentum 
fraction of the parton distribution function (PDF) value 
due to the variation of the jet masses between regions 
A to D. Furthermore, it is possible to test this assump- 
tion by using various Monte Carlo (MC) tools to extract 
-Rmass- We did so with four different tools. The results 
are summarized in Table [I] 

The impressive agreement between Sherpa and Mad- 
Graph when matching is employed leads us to use, in- 
stead of Eq. ([3]) , the estimate 



= 0.87. 



(4) 



The estimated number of background events within the 
data sample of 5.95 fb _1 is then 



QCD : 15 ±5, 
tt : 3 ± 1 . 



MC tools Matching 




-Rmass 


Sherpa 


Yes 


0.8 


8 ±0.03 


MadGraph 


Yes 


0.8 


6 ± 0.04 


MadGraph 


No 


0.76 ± 0.04 


Herwig 


No 


0.8 


6 ± 0.02 



TABLE I: The results for i? maaa [Eq. ©] from the differ- 
ent MC tools: Sherpa (1.2.3) [3(| with matching (jjjjjjjjj), 
MadGraph/MadEvent (4.4.56) [13] with MLM match- 
ing [H (jjjjjjjjj) to the Pythia package (2.1.4) Mad- 
Graph/MadEvent with no matching, and Herwig (6.520) [4CJ] 
with no matching We use the CTEQ6M PDF set [H| and 
FastJet (2.4.2) Q with anti-k t algorithm HH (AR = 1). 
Quoted errors are statistical only. 



The number of observed events was 32 [6j, which con- 
stitutes a deviation of 2.7er from the above expectation. 
Following the exercise performed in Ref. Q, the differ- 
ence between the 32 observed events and the mean value 
of Eq. ([5]) is translated to a cross section of 



NP 



10 ±4 fb, 



or, equivalently, 



AT — ~ NP /^-SM 

N b = a h /a b 



5±2. 



(6) 



(7) 



Below we obtain predictions from new physics scenarios 
for N b , which we will compare against Eq. ([?])• 

Additional data. Other top-related CDF and DO 
measurements, beyond A" and a b , do not show signifi- 
cant deviations from the SM predictions. (Interestingly, 
a recent DO measurement of the differential px distribu- 
tion of tt events hints towards some increase over the 
NLO SM prediction for p T ~ 300 GeV Q.) When we 
invoke new physics to account for the large value of A" , 
we will have to make sure that such new physics does not 
violate the constraints from other measurements. Specif- 
ically, we consider the following measurements: 

(i) The forward-backward tt production asymmetry for 
small invariant mass of the tt system [l| : 



J^tt _ Att 



A u {M tl < 450 GeV) = -0.116 ± 0.153 , (8) 



to be compared with the SM prediction 0], A" 
+0.040 ±0.006. 

(ii) The inclusive tt production cross section reported 
by the CDF Collaboration [45|,[46|: 



_ tt 

a i — inclusive 



= 7.50 ±0.48 pb, 



(9) 



This is to be 
a, = 7.2±0.4pb. 



(5) 



which is consistent with the DO result [47 
compared with the SM prediction 

We note that the results of Q| agree with other recent 
evaluations 0, IHoj . but are in some tension with [Hlj . 
We conservatively use this result, as that of [Hl[ would 
be less constraining given our framework. 
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(iii) The tt differential cross section, which for sim- 
plicity we choose to represent by the following large M t % 
bin [52J: 

a h = CT tf (700 GeV < M tf < 800 GeV) = 80 ± 37 fb , 

(10) 

to be compared with the SM prediction [2J, |51|, = 
80 ± 8 fb. The choice of this specific bin requires some 
explanation. 

• Since we focus on new physics which contributes to 
the tt cross section oc (M t i/ K) 2 relative to the SM, 
the corrections to lower M t j bins are less significant. 

• In the more recent study of [1], which was based on 
a larger sample, there is some discrepancy above 
800 GeV (note however that the data in [l| is not 
unfolded to the partonic level and so cannot be di- 
rectly used). Hence we choose to use the next-to- 
last bin given in (52| . 

In order to minimize the impact of NLO corrections 
to the new physics (NP) contributions, we normalize the 
new physics contribution to the SM one. We assume that 
the if -factors are universal, so that the NP/SM ratios at 
LO and NLO are the same. Since the highly virtual inter- 
mediate gluon in the SM process can be integrated out to 
give Oy , NP NLO contributions should be similar, up to 
small corrections of 0(a s ). Moreover, the parity invari- 
ance of QCD suggests that the same argument applies to 
G\ as well. 

Combining in quadrature the experimental and theo- 
retical uncertainties, we represent Eqs. (J9j) and (fT0|) as 
follows: 



Ni 



— unp /SMI <r n i 

AT — I NP /SMI < A r 

N h = \a h /a h <; 0.5 . 



(11) 



C e s for tt production. The basic assumption that 
we aim to test is that the source of the large value of A h * is 
new physics that is characterized by a mass scale A that 
is larger than M t i in all the measurements that we con- 
sider. (In particular, our Tevatron-related calculations 
are safe for A ^ 1 TeV.) In such a case, the new physics 
can be represented as a set of effective operators. These 
operators must lead from an initial uu state to a final tt 
state. (The contribution of dd — > tt at the Tevatron is 
at most 15% that of uu — > tt for M t f above 450 GeV, as 
relevant for the observables that we consider.) When ex- 
panding in inverse powers of the scale A, the leading NP 
contributions to top pair production appear at 0(1/A 2 ): 



IMI 



|M SM | 2 + 2Ke(M SM M^ P ) + 0(1/ A 4 



(12) 



Therefore, we should consider dimension-six operators 
that interfere with the SM amplitude. There are two 



such four-quark operators: 
1 



C 



off - ^2 ( C A®/ 

0\ = (uw 5 T a u)(i^ 5 T a t) , 
O v = (u^T a u)(t^T a t) . 



(13) 



Below, we consider the effects of these two operators on 
the forward-backward asymmetry and on the differential 
cross section in top pair production. We work only at 
leading order, using the MSTW PDF set [53j and run- 
ning of the strong coupling at leading order. We use fac- 
torization and renormalization scales given by the par- 
tonic center of mass energy. Note that all other possi- 
ble Lorentz structures (scalar, pseudoscalar, tensor and 
pseudotensor) and the other possible color contractions 
do not interfere with the SM amplitude. 

In addition to the four-quark operators, there is a chro- 
momagnetic dipole operator, 



C 



tg 

eff 



CtgV 

A 2 



{ia^T a t)G a ^ . 



(14) 



Here v is the vacuum expectation value of the Higgs field, 
reflecting the fact that the operator breaks SU(2). The 
corresponding chromoelectric dipole operator violates CP 
and, therefore, does not interfere with the SM amplitude. 
The interference of the chromomagnetic operator requires 
a chirality flip. Consequently, the corresponding operator 
involving the up quark is suppressed by m u and therefore 
negligible. Thus, among the dipole operators, Eq. ([T4")) is 
the only one that we need to consider. 

The interference of the Ct g term with the SM amplitude 
does not contribute to the forward-backward asymmetry. 
As concerns the contribution to the cross section, it falls 
like 1/Af t 2 j [32j]. We learn that while the Ct g term can 
affect the inclusive cross section, it does not affect A^, 
and its effect on Nh and Ni, is negligible. We therefore 
focus mainly on the effects of 0\ and 0\. See, however, 
additional discussion above Eq. (|2"3"1) . 

The forward-backward asymmetry. It is con- 
venient to represent the new physics effects on A u as 
follows: 



T NP 



(A") NP = 



V NP 



(15) 



where a± = <r(Ay > 0) ±a(Ay < 0) and Ay is the rapid- 
ity difference, Ay = y t — yi- Among the two operators of 
Eq. (p~5|) . only 0\ contributes to er_ . If this is the only 
NP operator, the NP contribution to A h l is 



tt\NP 



0.17 



A 2 



where A To v = A/TeV. Requiring that (^') NP 
0.1, we obtain 



4/A 



TcV 



2.4 ±0.7. 



(16) 
-0.4 ± 

(17) 
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Eq. (fl~7|) implies, in turn, 
(Aff p ~ +0.10 ±0.03 



A? 



-0.14 ±0.04, (18) 



about 1.7 'a higher than the experimental result in Eq. (jH). 
In addition, Eq. ^7} predicts (A") NP - +0.21 ± 0.06, 
1.5(7 too large and (A**(Ay > 1)) NP - +0.55 ± 0.15, 
within one standard deviation from the measurement [lj . 
On the other hand, the 0\ operator does not affect the 
cross section at 0(1/A 2 ) and, in particular, cannot en- 
hance the boosted jets cross section. The contribution of 
the next order in 1/A 2 to the forward-backward asym- 
metry is subdominant, and, using the one sigma lower 
bound of Eq. (jTTJ) , saturates the constraint from Nh in 
Eq. (HU). 

Eq. (|17[) provides an upper bound on the scale of new 
physics. We use naive dimensional analysis (NDA) to 
derive an upper bound on c 1 



■Ai 



(19) 



Combining this upper bound with the one sigma lower 
bound in Eq. (|17l). we obtain 



A < 10 TeV . 



(20) 



The upper bound on A in Eq. (|20|) implies that, if 
a heavy axigluon is to provide a pcrturbative explana- 
tion to the large asymmetry in Eq. ([1]) , then new physics 
effects should be observed early on at the LHC. In par- 
ticular, given that the LHC will directly explore energy 
scales close to A , then the tt production cross section 
should be significantly enhanced at high M t i (see [13] for 
more details). 

To substantiate this statement, we perform the follow- 
ing exercise. We note that the 0\ operator does modify 
the cross section at 0(1/A 4 ) via the |Mnp| 2 term. We 
plot in Fig. [1] the differential tt cross section as a func- 
tion of M t i a t the LHC for the case where the SM is 
augmented by only the 0\ operator, with the coupling 
of Eq. (fl~7| (the distribution at the Tevatron is also de- 
picted for comparison). Of course, at this order there are 
many more operators that affect the cross section, either 
via |M NP | 2 for C(l/A 2 ) operators , or via T^(Msm-^np) 
for C(l/A 4 ) operators. In Ref. [54] it is shown, how- 
ever, that there can be no fine-tuned cancellations be- 
tween these other contributions and the one that we con- 
sider. Thus our calculation illustrates the size of the ef- 
fects that should be expected at the LHC. We learn that 
at M t i ~ 1.5 TeV, we should expect an enhancement by 
a factor ~ 5 compared to the SM. When applied to the 
Tevatron, the same exercise gives an enhanced boosted 
jets production cross section of JVj ~ 2, which is 1.5er 
from the mean value in Eq. ([7]). Fig. [2] depicts the result- 
ing pt distribution at the Tevatron. 

We stress that the recent measurement of the differ- 
ential tt forward-backward asymmetry predicts a more 



pronounced deviation from the SM at the LHC than the 
previous inclusive asymmetry measurement. This is illus- 
trated in Fig. [1] by the difference between the solid and 
dashed curves (and the respective shaded regions) and 
the dashed-dotted curve representing the SM prediction. 
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FIG. 1: The differential cross section of top pair production 
as a function of M t j at the LHC at 7 TeV (top) and the 
Tevatron (bottom), calculated at leading order. The dashed- 
dotted curve corresponds to <r£o • The solid curve corresponds 
to a SM+ ° A , where the new physics coupling is set by the cen- 
tral value of A™. The dashed curve corresponds to <j SM+ ° A , 
where the new physics coupling is set by the central value of 
the inclusive asymmetry A u . The shaded regions around the 
two upper curves depict the one sigma ranges of the corre- 
sponding measurements. 

The tt production cross section. Among the 



two operators of Eq. 1(13)1. only Oy contributes to the 
inclusive cross section (ov), to the cross section at large 
M t i (o'h) and to the production cross section of boosted 
jets (o-fe) at 0(1/A 2 ): 



Ni 
N h 



0.24 c s v /A 2 TeV , 
0.76 4/A 2 oV , 

ToV ■ 



1.5 c^/A 



(21) 



This equation, where relevant, is consistent with previous 
results in the literature ([12] and references therein). In- 
dependently of the value of the coupling, our framework 
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predicts 



(22) 



This ordering of the size of the effects reflects the fact that 
each of these three measurements samples a different M t i 
region; the closer this region is to A, the larger the effect. 

The relation between N b and Ni and between Nh and 
Ni can be modified by the presence of the chromomag- 
netic dipole operator in Eq. (fl4|) . if |c tg | is not much 
smaller than \c v . However, since the Ct g term does not 
affect the cross section at high invariant mass M% 3> m 2 , 
the relation between Nf, and Nh is insensitive to it. The 
bound on Nh in Eq. ([TTj) then leads to an upper bound 
on the enhancement of boosted jets production: 



Obs. 


Def. 


Experiment 


Standard Model 


New physics 


At 


Eq. © 


+0.475 ± 0.114 


+0.09 ± 0.01 


Input 


Af 


Eq. © 


-0.116 ±0.153 


+0.040 ± 0.006 


+0.16 + 0.04 


CTi 


Eq. © 


7.50 ± 0.48 pb 


7.2 + 0.4 pb 


Input 


Oh 


Eq. (TO]) 


80 ± 37 fb 


80 + 8 fb 


Input 


O b 


Eqs. (1214)) 


12 ±4 fb 


2.0 + 0.2 fb 


< 3.2 fb 



TABLE II: Effects from new physics of C(l/A 2 ) on top- 
related observables. The first column gives the list of observ- 
ables, and the second the equation where they are defined. 
We use At, &i and on to fix, or constrain, the new physics 
parameters. The experimental value quoted for at is based 
on our theoretical interpretation of the data. 



N b < 0.8 . 



(23) 



well below our estimate of Eq. CZ]). We conclude that one 
of the following must hold: 

• The estimate of Eq. CB is wrong because of either 
experimental or QCD effects. 

• New physics explains Eq. (O, but it is characterized 
by a scale that is < 1 TeV. 

• Heavy new physics explains Eq. CB, but 0(1/ A 4 ) 
terms play an important role [541 ] . 

• The reported excess in events with two boosted 
massive jets does not originate from top quarks. 

Conclusions. 



tt forward-backward asymmetry at large M t i , , 



The recent CDF measurement of the 

At, shows 

a deviation higher than 3<r from the SM prediction. The 
recent CDF measurement of ultra-massive boosted jets, 
<7b, shows a deviation of order 2.7c from a SM calculation 
augmented with an estimate of QCD background based 




0.001 
350 



550 



FIG. 2: The differential cross section of top pair production as 
a function of pr at the Tevatron, calculated at leading order. 
The color and curve conventions are the same as in Fig. [1] 
The vertical line corresponds to the lower pr cut used in the 
analysis of 



on data and on several simplifying assumptions that we 
test with various MC tools. 

We investigated whether these effects can be accounted 
for within a large class of new physics models. This class 
of models is defined by a mass scale above the scales di- 
rectly explored by these CDF measurements, and a dom- 
inant effect coming from interference between the Stan- 
dard Model and new physics amplitudes. 

Within this framework, we find that there is a sin- 
gle four quark operator that contributes to the asym- 
metry, which is the axial vector, color octet, operator 
Q\ = (u~f^T a u)(i^^T a t). There is a different single 
four quark operator that modifies the differential cross 
section at high tt invariant mass, which is the vector, 



color octet, operator Of 



(u-y IJ T a u)(t 1 ^T a t). This 



means in particular that there is no model independent 
relation between the forward-backward asymmetry and 
the boosted jets cross section. Note that we focus on 
these operators at tree level, and so do not discuss their 
contribution to dijet production at the LHC 



12] 



Our numerical results are summarized in Table HU If 
0\ accounts for the high value of At , then the asymme- 
try at low invariant mass is about 1.7a high compared 
to the CDF measurement. One should expect a striking 
enhancement of tt production at high M t i at the LHC. 

If Oy is to be consistent with constraints from the in- 
clusive and differential cross sections, then it cannot en- 
hance the boosted tops cross section by more than 60%. 
Furthermore, O v is restricted to be significantly smaller 
than the contribution of 0\ implied by the tt asymme- 
try. This means that a chiral model cannot consistently 
reproduce the asymmetry. 

The above conclusions are related to the fact that the 
interference effects of heavy new physics with the SM 
scale roughly as (M t f/A) 2 relative to the SM. Conse- 
quently, they do not differentiate between the low and 
high M t i regions enough to avoid tension with the data. 

The conclusion concerning the ultra-massive boosted 
tops is that 0(1/A 2 ) effects do not explain the discrep- 
ancy of the data with our theoretical estimate of the SM 
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contribution. Perhaps the explanation does not involve 
new physics: The deviation is below 3tr and might dis- 
appear with better experimental accuracy, or it could 
be that QCD effects that are unaccounted for in the 
various MC tools play a role. If the deviation is re- 
lated to new physics, then either the new physics is be- 
low the TeV scale and cannot be represented by effec- 
tive higher-dimension operators, or the contribution of 
|Mnp| 2 oc 1/A 4 is significant, bringing into the analysis a 
richer set of operators and a sharper distinction between 
the low and high M t i regions. 

The LHC will explore ti production at higher energy 
scales. Whether the scale A is within its direct reach 
or just beyond it, new physics effects are expected to 
be large. The Tevatron, on the other hand, has better 
access to the qq — > ti process which, via observables such 
as the forward-backward asymmetry, can close in on the 
detailed structure of new physics. The combination of 
Tevatron and LHC measurements is likely to shed light 
on the top-related puzzles very soon. 

Acknowledgments. We thank Raz Alon, Yochay 
Eshcl, Alex Kagan, Fabrizio Margaroli, Pekka Sinervo 
and George Sterman for useful discussions. YN is the 
Amos de-Shalit chair of theoretical physics and supported 
by the Israel Science Foundation (grant #377/07), the 
German-Israeli foundation for scientific research and de- 
velopment (GIF), and the United States- Israel Binational 
Science Foundation (BSF) . GP is the Shlomo and Michla 
Tomarin career development chair and is supported by 
the Israel Science Foundation (grant #1087/09), EU-FP7 
Marie Curie, IRG fellowship, Minerva and G.I.F., the 
German-Israeli Foundations, and the Peter & Patricia 
Gruber Award. 



* Electronic 
' Electronic 

* Electronic 
§ Electronic 
" Electronic 
** Electronic 
*' Electronic 

Electronic 



address 
address 
address 
address 
address 
address 
address 
address 



kfir . blum@weizmann.ac.il 
cedric.delaunay@weizmann.ac.il 
oram.gedalia@weizmann.ac.il 
yonit . hochberg@weizmann.ac.il 
seung.lee@weizmann.ac.il 
yosef.nir@weizmann.ac.il 
gilad.perez@wcizmann.ac.il 
yotam.soreq@weizmann.ac.il 



[1] T. Aaltonen et al. [CDF Collaboration], Phys. Rev. D 

83, 112003 (2011) [arXiv:1101.0034l [hep-ex]]. 
[2] L. G. Almeida, G. F. Sterman and W. Vogelsang, Phys. 

Rev. D 78, 014008 (2008) [arXiv: 0805 . 18851 [hep-ph] ] . 
[3] M. T. Bowen, S. D. E llis and D. Rainwater Phys. Rev. 

D 73, 014008 (2006) [arXiv:hep-ph/0509267j . 
[4] O. Antunano, J. H. Kuhn and G. Rodrigo, Phys. Rev. D 

77, 014003 (2008) [arXiv:0709. 16521 [hep-ph]]. 
[5] CDF Collaboration, CDF note 10199, January 18, 2011. 
[6] CDF Collaboration, CDF note 10234, January 5, 2011. 
[7] Y. Eshel, O. Gedalia, G. Perez and Y. Sor eq, Phys. Rev. 

D 84, 011505(R) (2011) [arXiv: 110 1.28981 [hep-ph]]. 
[8] N. Kidonakis and R. Vogt, Phys. Rev. D 68, 114014 



[9 
[10 

[11 

[12 

[is; 
[14; 

[is; 
[is; 
[17; 
[is; 

[19 
[20 
[21 
[22 
[23 
[24 
[25 
[26 

[27 
[28 
[29 
[30 
[31 
[32 

[33 
[34 
[35 
[36 

[37 
[38 



(2003) [arXiv:h ep-ph/0308222] . 

K. Che ung and T. C. Y uaiTTPhys. Rev. D 83, 074006 
(2011) |arXiv:1101. 14451 [hep-ph]]. 

C. Delaunay, O. Gedalia, S. J. Lee, G. Perez and E. Pon- 
ton, larXiv:1101.2902l [hep-ph] . 

J. Cao, L. Wang, L. Wu and J. M. Yang. larXiv:1101.4456l 
[hep-ph]. 

Y. Bai, J. L. Hew ett, J. K aplan and T. G. Rizzo, JHEP 
1103, 003 (2011) [arXivillOl. 52031 [hep-ph]]. 
J. Shel ton and K. M. Z urek, Phys. Rev. D 83, 091701 
(2011) |arXiv:1101. 53921 [hep-ph]]. 

E. L. Berger, Q. H. Cao, C. R. Chen, C. S. Li 
and H. Zhang, Phys. Rev. Lett. 106, 201801 (2011) 
[arXivillOl. 56251 [hep-ph]]. 

B. Bhattacherjee, S. S. Biswal and D. Ghosh, Phys. Rev. 

D 83, 091501 (2011) |arXiv:1102.0545l [hep-ph]]. 

P. Ferr ario and G. Rod rigo, Phys. Rev. D 78, 094018 

(2008) |arXiv:0809.3354l [hep-ph]]. 

P. Ferr ario and G. Rod rigo, Phys. Rev. D 80, 051701 

(2009) iarXiv:0906.554T1 [hep-ph] j. 

S. Jung, H. Murayama, A. Pierce and J. D. Wells, Phys. 

Rev. D 81, 015004 (2010) |arXiv:0907.4TT2l [hep-ph]]. 

K. Cheung, W. Y. Keung and T. C . Yuan, Phys. Lett. B 

682, 287 (2009) [arXiv: 0908. 25891 [hep-ph]] . 

P. H. Fram pton, J. Shu and K . Wang, Phys. Lett. B 683, 

294 (2010) [arXiv:0911. 29551 [hep-ph]]. 

J. Shu, T. M. P. Tait and K. Wang, Phys. Rev. D 81, 

034012 (2010) |arXiv:0911.3237l [hep-ph]]. 

A. Arhrib, R. Benbrik and C. H. Chen, Phys. Rev. D 82, 
034034 (2010) |arXiv:0911.4875l [hep-ph]]. 

I. Dorsner, S. Fajfer, J. F. Kamenik and N. Kosnik, Phys. 

Rev. D 81, 055009 (2010) |arXiv:0912.0972l [hep-ph]]. 

J. Cao, Z. Heng, L. Wu and J. M. Yang, Phys. Rev. D 

81, 014016 (2010) [arXiv:0912.1447l [hep-ph]]. 

V. Barger, W. Y. Keung and C. T. Yu, Phys. Rev. D 81, 

113009 (2010) arXiv:1002. 10481 [hep-ph]]. 

Q. H. Cao, D. McKeen, J. L. Rosner, G. Shaughnessy 

and C. E. M . Wagner, Phys. Rev. D 81, 114004 (2010) 

[arXiv:1003.346T1 [hep-ph]]. 

B. Xiao, Y. k. Wang and S. h. Zhu, Phys. Rev. D 82, 
034026 (2010) arXiv:1006. 25101 [hep-ph]]. 

R. S. Chivukula, E. H. Simmons and C. P. Yuan, Phys. 
Rev. D 82, 094009 (2010) |arXiv: 1007.02601 [hep-ph]]. 

C. H. Chen, G. Cvetic and C . S. Kim, Phys. Lett. B 694, 
393 (2011) [arXiv:1009.4"T65l [hep-ph]]. 

E. Alvarez, L. Da Rold and A. Szynkman, JHEP 1105, 
070 (2011) [arXiv:1011.6"557l [hep-ph]]. 

D. W. Jung, P. Ko, J. S. Lee and S. h. Nam, Phys. Lett. 
B 691, 238 (2010) arXiv:0912.TT05l [hep-ph]] . 

C. Degrande, J. M. Gerard, C. Grojean, F. Maltoni and 
G. Servant, JHEP 1103, 125 (2011) |arXiv: 10 10.63041 
[hep-ph]]. 

D. W. Jung, P. Ko and J . S. Lee, Phys. Lett. B 701, 248 
(2011) |arXiv:1011. 59761 [hep-ph]]. 

V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 

100, 142002 (2008) (arXiv:0712.0851l [hep-ex]]. 

T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett. 

101, 202001 (2008) |arXiv:0806. 24721 [hep-ex]]. 

T. Gleisberg, S. Hoeche, F. Krauss, M. Schonherr, 
S. Schu mann, F. Siegert and J. Winter, JHEP 0902, 007 

(2009) [arXiv:0811.4622l [hep-ph]]. 

J. Alwall et al, JHEP 0709, 028 (2007) |arXiv:0706. 23341 
[hep-ph]]. 

S. Hoeche, F. Krauss, N. Lavesson, L. Lonnblad, 



7 



M. Mangano, A. Schalicke and S. Schumann, [48] 

[arXiv:hep-ph/ 0602031 . 
[39] T. Sjostran d, S. Mrenna and P. Z Skands, JHEP 0605, [49] 

026 (2006) [arXiv:h ep-ph/06031 75]. 
[40] G. Corcella et al, [arXiv:hep-ph/0210213] . 
[41] J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, [50] 

P. M. Nadolsky and W K. Tung, JHEP 0207, 012 (2002) 

[arXiv:hep-ph/0201195] . [51] 
[42] M. Ca cciari and G. P. Salam , Phys. Lett. B 641, 57 

(2006) [arXiv:hep-ph/0512210| , |http : //f ast j et . f r/| 
[43] M. Cac ciari, G. P. Salam and G. Soyez, JHEP 0804, 063 [52] 

(2008) [arXiv:0802. UMI [hep-ph]], 
[44] V. M. Abazov et al. [DO Collabo ration], Phys. Lett. B [53] 

693. 515 (2010) arXiv: 1001. 1900 1 [hep-ex]]. 
[45] CDF Collaboration, CDF note 9913, October 19, 2009. 
[46] C. Schwanenberger, for the CDF Collaboration and [54] 

the DO Collaboration, larXiv:1012.23T9l [hep-ex] . 
[47] DO Collaboration, DO note 6037, May 28, 2010. 



D 82, 114030 (2010) 



Mangano, P. Nason and 
(2008) |arXiv: 0804.2800 



N. Kidonakis, Phys. Rev. 
[arXiv: 1009.49351 [hep-ph]]. 
M. Cacciari, S. Frixione, M. L. 
G. Ridolfi, JHEP 0809, 127 
[hep-ph]]. 

U. Langenfeld, S. Moch and P. Uwer, Phys. Rev. D 80, 

054009 (2009) |arXiv:0906.5"273l [hep-ph]]. 

V. Ahrens, A. Ferroglia, M. Neubert, B. D. Pecjak and 

L. L. Yang, JHEP 1009, 097 (2010) |arXiv: 1003.58271 

[hep-ph]]. 

T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett. 
102, 222003 (2009) )arXiv:0903.2850l [hep-ex]]. 
A. D. Martin, W. J. Stirling, R. S . Thome and G. W att, 
Eur. Phys. J. C 63, 189 (2009) |arXiv:0901.0002l [hep- 
ph]]. 

C. Delaunay, O. Gedalia, Y. Hochberg, G. Perez and 
Y. Soreq, larXiv:1103.2297l [hep-ph] . 



